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TECHNICAL NOTE NO. 1792 


STUDY BY THE IRANDTL-GLAjQERT METHOD OF COMPRESSIBILITY 
EFFECTS AND CRITICAL MACH NUMBER FOR ELLIP SOIDS OF 
VARIOUS ASPECT RATIOS AND THICKNESS RATIOS 
By Robert V. Hess and Clifford S. Gardner 

SUMMARY 


By the use of a form of the Prandtl-Glauert method that is valid 
for three-dimensional flow problems, the value of the maximum incremental 
velocity for compressible flow about thin ellipsoids at zero angle of 
attack is calculated as a function of the Mach number for various aspect 
ratios and thickness ratios. The critical Mach numbers (within the 
accuracy of the Prandtl-Glauert method) of the various ellipsoids are 
also determined. The results indicate an increase in critical Mach 
number with decrease in aspect ratio which is large enough to explain 
experimental results on low-aspect -ratio wings at zero lift. 


INTRODUCTION 


Recent tests (references 1 and 2) have shown that an appreciable 
Increase in the critical Mach number, together with other improvements 
of the aerodynamic characteristics at supercritical Mach numbers, results 
from the use of wings of very low aspect ratio- These improved charac- 
teristics have been somewhat qualitatively ascribed to "three-dimensional 
relief, although no quantitative theoretical discussion has yet been 
provided- 

In the present paper an effort is made to provide such a study by 
considering the flow, at zero angle of attack, about a series of thin 
ellipsoids of various aspect ratios and thickness ratios. Ellipsoids 
were chosen because they are amenable to calculation. Although they 
differ appreciably from the wings of reference 1, which had an NACA 
0012 airfoil section and rectangular plan form, ellipsoids should never- 
theless show similar aspect-ratio effects- The calculations were made 
for ellipsoids of thickness ratios 0.10, 0-15, and 0-20, and for the 
entire range of aspect ratios from the elliptic cylinder to the 
ellipsoid of revolution. 

The compressibility effects were computed by the use of a form of 
the Prandtl-Glauert method that Is valid for three-dimensional flow 
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problems . The method, has been given "by GSthert (reference 3) without, 
however, very explicit mathematical proof* Another correct statement 
of the three-dimensional form of the Prandtl-Glauert method was given 
earlier "by A* Busemann in reference k-, where, however, no formulas were 
given* Since the methods that have been commonly used (see, for example, 
references 5, 6, arid 7) are applicable only to two-dimensional problems, 
a detailed proof of the method correct for three-dimensional flow is 
included in the appendix* A brief discussion of the accuracy of the 
Prandtl-Glauert method, as applied to ellipsoids, is. also given. This 
work was completed in April 19^6* 

SIMBOIS 
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free-stream velocity 

velocity of sound in free stream 

free-stream Mach number (u/c) 

ratio of specific heats (7 = lA for air) 


rectangular coordinates 

thin body 

velocity potential 

x-, y-, and z-ccmponents of incremental 
velocity for compressible flow about B 

body obt ain ed by stretching B in direction 
of x-axis by the factor l/p 

x-, y-, and z-ccmponents of incremental 

velocity for incompressible flow about B 1 

maxInriTm semi chord of ellipsoid 



semispan of ellipsoid 

maximum semithickness of ellipsoid 
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A 

aspect ratio ( A = ^ = 4 A 
y «ab it ay 

- %ax 
u ” U 


u(M) 

value of u when the Mach number is equal to M 

u(0) 

value of u for incompressible flow (M = 0) 

6 

/ Thickness'N 

thickness ratio 

y Chord J 


value of ratio of incremental velocity to 
free-stream velocity for compressible 
flow having Mach number M about a body 
having thickness ratio e 


value of ratio of incremental velocity to free- 
stream velocity for incompres Bible flow about 
a body having thickness ratio e 


Subscript: 

max maximum value 


METHODS OF CALCULATION 


The Pr andtl -Glauert method for three-dimensional flow . - The 
Prandtl-Glauert method is used in the present paper in the following 
form: 


The incremental velocities at a point P of a three-dimensional 
compressible flow field about a thin body B may be obtained in t hr ee 
steps: 

(l) The x-coordinates of all points of B are increased by 
the factor l/(3, where 

P = \]l - M 2 

and where the x-axis is in the stream direction. This transformation 
takes B into a stretched body B 1 . 
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(2) The inc remental velocities u 1 , v 1 , w', in the direction of 
the x-, y-, and z-axes, respectively, at the point P 1 in the flow 
field of B 1 corresponding to the point P in the flow field of B 
are calculated as though B 1 were in an incompressible flaw having the 
same free-stream velocity as the original compressible flow- 

( 3 ) The values u, v, and w of the incremental velocities at the 
point P in the compressible flow field of the original unstretched 
body are then found by the equations 


1 , 
u = — u’ 

3 2 

1 , 

V = — V 

3 


w = — w' 
3 


A derivation of this form of the Prandtl-Glauert method is given 
in the appendix. The method in essentially this form has been given by 
Gothert (reference 3 ) without, however, a very clear proof. (GcSthert 
prefers to shrink the lateral coordinates of the body by the factor 3 
rather than to expand the coordinate in the stream direction by the 
factor l/ 3 j obviously the two procedures lead to the same result.) 

Prandtl (reference 5) and Ton Kerman (reference 6 ) state the method in 
a form that is valid for two-dimensional flows but in general is incorrect 
for three-dimensional flows. Goldstein and Young (reference 7) also give 
a discussion leading to results that are correct only for two dimensions. 

A discussion of the reasons for the failure of these commonly used 
methods for three-dimensional flow problems is included in the appendix* 

Calculation of Incremental velocity for compressible flow about 
ellipsoids . - In order to dete rmine , by the Prandtl-Glauert method, the 
incremental velocity on the surface of an ellipsoid having semiaxes a, 
b, and c (where a • is the length of the semiaxis in the stream direc- 
tion), the incremental velocity is calculated for a stretched ellipsoid 

hav ing semiaxes a', b, and c (where a' = ^ in an incompressible 

‘3 

flow having the same stream velocity) and the result is multiplied 

by l/3 2 . For incompressible flow about the stretched ellipsoid, the 
velocity potential on the surface of the ellipsoid is given by 
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where 


dr 



6X 

(a' 2 + x) ^(a ' 2 + x) (b 2 + x) ^c 2 + x) 


(See, for example, reference 8.) The incremental velocity at x = 0 
(half-chord, line on the stretched, ellipsoid in incompressible flow) is 
then given by 


u 


1 



U 


This value is the maximum value of u 1 (reference 8) and evidently 
is the same at a ll points on the half-chord l in e. The incremental 
velocity at the half -chord line for the compressible flow about the 
original unstretched ellipsoid is given by 


u 



1 a o 

p 2 2 - a Q 


TJ 


( 1 ) 


Various formulas are necessary for the evaluation of the integral a Q 
when a 1 > b > c, b > a' > c, or a 1 > b = c (ellipsoid of revolution). 

Far a 1 > b > c, the value of d Q is given by the formula 


2a'bc 

% — 

(a' 2 - b 2 ) \/a' 2 - c 2 

where F and E are incomplete elliptic integrals of the first and 
second kind, respectively, defined as follows: 


(F - E) 


( 2 ) 


19 


F = 


dt 


Vl - k 2 sin 2 i|r 


E = 


f \jl - k 2 sin 2 V di{f 
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where 


/ a' 2 - b 2 
^ \j a' 2 - c 2 

and 


Va' 2 ~ c 2 

sin cp = i 

a 


For b > a 1 > c, the value of a Q Is given by the formula 


% = 


2a'bc ^b 2 - c 2 


G > 8 - °' 2 ) G ' 2 - c 2 ) L 


,2 2 ' 

a - c 

1 -1 I F 


-J2 2 

^b - c 


2o £ 


.2 


(3) 


where F and E are defined by the foregoing formulas but 


k = 


f t 2 - a* 2 
lb 2 - c 2 


add 


sin cp = 




Equation (2) is derived from the first equation given in equa- 
tions ( 5 * 13 ) of reference 8 by substituting a 1 for a and by using 
the expression for k in terms of a 1 , b, and c. Equation ( 3 ) ia 
derived from the second equation given in equations ( 5 * 13 ) of reference 
by interchanging a and b, substituting a 1 for a , and using the 
expression for k in terms of a 1 , b, and c. 

For a 1 > b = c (ellipsoid of revolution), Oq is given by the 
equation 


a D = a'b 2 


dX 


l 0 (a'2 + X 2 )3/2(b 2 + x) 
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which resolves into 


where 




If this value for Oq is substituted in equation (l) , the incremental 
velocity at the half -chord line for the ellipsoid of revolution is found 
to he 


u 


1_ 

P 2 


, 1 + e 

log e " 2e 

2e " 1 + e 

“ logQ 

1 - e d 1 - e 


U 


The limiting case of infinite aspect ratio (elliptic cylinder) was 
treated by the use of formulas for the ellipse in two-dimensional flow 
(reference 9) • 

Calculation of the critical Mach number . - For flow about a two- 
dimensional body , the free -stream Mach number for which sonic speed is 
first reached at seme point on the surface is called the critical Mach 
number , because of the development of a shock and the accompanying 
deterioration of the aerodynamic characteristics shortly after this 
Mach number is exceeded. If an extension of the definition of the 
critical Mach number to the general three-dimensional body is desired, 
the definition appears, at first sight, to require same revision, since 
for the general three-dimensional case the shock formation on parts of 
the body may occur along lines yawed with respect to the free-stream 
velocity (reference 10) • The boundary lines of the supersonic regions 
(sonic lines) must, however, always contain a portion that is normal 
to the free-stream velocity and thus the definition of the critical 
Mach number for the two-d i mensional case (infinite unyawed cylinder) 
may be extended without modification to the general three-dimensional 
case- For the special three-dimensional case of the infinite yawed 
cylinder, the portions of 'the boundary lines of the supersonic regions 
that are normal to the free-stream velocity are represented by two 
points at infi n ity, only one of which (at the downstream end) has the 
necessary qualities for accumulating disturbances, that is, for shock 
formation. It can be seen that for the general three-dimensional body 
the critical Mach number, although still defined in the same mann er as 
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for the two-dimensional body, largely loses its critical significance, 
since the shock may begin, to form over only a very sma3 1 part of the 
surface, so that its occurrence does not necessarily imply an imminent 
deterioration of the aerodynamic characteristics of the body. 

For the special case of the unyawed ellipsoid considered in the 
present paper, however, no appreciable analysis of shock formation or 
shock extent along the lines Just Indicated seems to be required. As 
is shown in the section "Calculation of the Incremental Velocity for 
Compressible Flow about Ellipsoids", the maximum velocity for an unyawed 
ellipsoid is in the stream direction and occurs simultaneously at all 
points along the half-chord line. Sonic velocity is thus reached 
simultaneously along a line that extends across the entire span of the 
body and is normal to the stream direction. These conditions also 
exist in the case of the unyawed infinite cylinder, that is, the two- 
dimensional body. 

The critical Mach number of the ellipsoid, within the accuracy of 
the Prandtl-Glauert method, was accordingly determined by solving 
graphically the equation 


7 - 1 o 

1 + ~2~ if 

1 

7 + 1 

2 

where u(m) is the ratio of the incremental velocity at the half-chord 
line to the stream velocity at the Mach number M. 

Accuracy of the Prandtl-Glauert method . - The Prandtl-Glauert method 
is based on the assumption of small perturbations. Consequently, near 
the nose of the ell i psoids discussed in the present paper, where the 
assumption of small perturbations is violated, the results given by the 
Prandtl-Glauert method cannot be expected to be reliable- More reliable 
values, however, should be obtained for the maximum incremental velocity, 
which occurs at the half -chord line. The accuracy of the Prandtl-Glauert 
approximation far the maximum incremental velocity may be estimated by 
comparison with more exact solutions of the compressible flow problem. 

An iteration method in which the Prandtl-Glauert method is used as the 
first approximation has been proposed by Busemann (reference ll) • The 
first and second approximations have been calculated by Hantzsche and 
Wendt for the elliptic cylinder (reference 12) and by Schmieden and 
Kawalki for the ellipsoid of revolution (reference 13) * Calculation 
of the maximum incremental velocity for the elliptic cylinder having 
thickness ratio 0.20 by a formula for the second approximation given 
in reference 11 shows that the value given by the Prandtl-Glauert method 
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at a Macli number of 0*8 is almost 20 percent lower than tbe value given 
by the second approximation- For the ellipsoid of revolution, however, 
the value of the max i mum incremental velocity given by the Prandtl-Glauert 
method agreed with the value given by the second approx im ation to within 
5 percent at a Mach number of 0-8 for thickness ratios up to 0-30- Although 
the second approximation is not the exact solution, it indicates that the 
error involved in using 'the Prandtl-Glauert method to estimate the maximum 
incremental velocity for ellipsoids having a given thickness ratio is 
greatest for the limiting case of the elliptic cylinder (A = 00 ) and 
very small for the ellipsoid of revolution, which has a very low aspect 
ratio- The error may be expected to be intermediate in magnitude for 0 
intermediate values of the aspect ratio and to decrease with aspect ratio. 
The reduction of error of the Prandtl-Glauert method with a decrease in 
aspect ratio was to be expected, as the incremental velocities are smaller 
for ellipsoids having low aspect ratio. 


RESULTS MD DISCUSSION 


Results . - Figures 1 to 3 show the value of the velocity 


ratio u = at the half -chord line plotted against the Mach number 

for ellipsoids at zero angle of attack for various aspect ratios and sec- 
tion thickness ratios equal to 0.10, 0-15# and 0-20- In the same figures 
the sonic velocity boundary having the equation 


u 



- 1 


is plotted for air. The abscissa of the intersection of this boundary line 
with the curve of u plotted against M for any aspect ratio is the 
critical Mach number (within the accuracy of the Prandtl-Glauert method) • 

In order to show the effect of compressibility more directly, the 

ratio of maximum incremental velocity for compressible flow to 

the maximum, incremental velocity for incompressible flow for the same 
free-stream velocity is plotted against the Mach number in figures 4 
to 6 for the same aspect ratios and thickness ratios- Similar curves 
for the ellipsoid of revolution, which is a special case of the ellipsoid 
having three unequal axes, are plotted for the same thickness ratios 
in figures 1 to 6- Figure 7 presents curves of critical Mach number 
against aspect ratio for thickness ratios of 0-10, 0-15, and 0-20- 

Three-dimensional relief . - It may be seen from figures 1 to 3 that 
the three-dimensional relief, that Is, the difference between the velocity 
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on the ellipsoid and the velocity on the corresponding ellipsoid of infinite 
aspect ratio (elliptic cylinder) , Increases with a decrease In the aspect 
ratio. This increase has two causes: 

(1) For a flow at amal 1 values of M (incompressible flow), the 
relief effect increases with a decrease in the aspect ratio. 

( 2 ) For larger values of M (compressible flow), an additional 
relief effect occurs with a decrease in the aspect -ratio because of the 
fact that the compressibility effect (increase of incremental velocity 
with an increase in the Mach number) decreases with a decrease in the 
aspect ratio. (See figs. 4 to 6 .) It may be seen that this additional 
three-dimensional relief increases most rapidly at high Mach numbers. 

Figures 1 to 6 show that the compressibility effect on the maximum 
incremental velocity is greatest far A equal to infinity (infinite 
elliptic cylinder) and is smallest for the ellipsoid of revolution. The 
compressibility effect on the maximum, incremental velocity far the elliptic 

cylinder is proportional to -, which is in agreement with the 

Vl - 

usual form of the Prandtl-Glauert method in two dimensions. The com- 
pressibility effect on the maximum incremental velocity for the ellipsoid 
of revolution is small, in comparison with that of the elliptic cylinder. 

In fact, as the thickness ratio of any type of body of revolution ' 
approaches zero, the compressibility correction factor approaches unity, 
for in this limit the incremental velocity in incompressible flow is 
proportional to the square of the thickness ratio, so that the effect of 
stretching the body (first step of Prandtl-Glauert method, see appendix) 
is exactly compensated for by the multiplication of the incremental 

velocities by l/p 2 (third step of the Prandtl-Glauert method). For 
ellipsoids of practical thickness ratios, however, the incremental 
velocity varies more slowly than the square of the thickness ratio. 

The compressibility effect for the ellipsoid of revolution (figs- 4 to 6 ) 
is thus appreciable at high Mach numbers. For example, for a thickness 
ratio of 0-20 and at a Mach number of 0.8, the compressibility effect 
amounts to about 30 percent of the incremental velocity in incompressible 
flow. 


The effect of the thickness ratio on the three-dimensional relief 
may be seen by a comparison of figures 1, 2, and 3* From figure 1 it 
may be seen that, for a thickness ratio of 0.10, at a Mach number of 0-75, 

the maximum incremental velocity for A = 2 is 76 percent of the maximum 

incremental velocity for A = 00 . From figure 3 ) on the other hand. It 
may be seen that, for a thickness ratio of 0.20, at a Mach number of 0-75, 

the maximum incremental velocity for A = 2 is 75 percent of the maximum 

incremental velocity for A = co. Thus, an increase in the thickness ratio 
causes only a very emal 1 increase in the three- dimens ional relief. 
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Critical Mach, number . Figures 1, 2, 3 * and 7 indicate that an increase 
in the critical Mach number of an ellipsoid at zero lift may he obtained 
by decreasing the aspect ratio. For example, for- ellipsoids having a 
thickness ratio of 0*10, a decrease in the aspect yatio from 00 to 2 
causes the critical Mach number to increase from O .827 to 0 . 857 (a Mach 
number increase of 0.Q3). For a thickness ratio of 0-20, a decrease in 
the aspect ratio from ® to 2 causes the critical Mach number to increase 
from 0.7^-l to O .783 (about 0*04). Although e lli psoids having greater 
thickness ratio have lower critical Mach n umb ers, a decrease in the 
aspect ratio is slightly more effective in inc re asing the critical Mach 
numbers for ellipsoids of greater thickness ratio. Figure 7 indicates 
that only a large reduction in aspect ratio will cause a significant rise 
of the critical Mach number. 

Comparison with test results on low aspect ratio wirigs .- Figure 6 
of reference 1 shows the minimum drag coefficient ^Cp for zero lift) 

plotted against the Mach number for wings having an NACA 0012 section 
and various aspect ratios. The critical Mach number for any aspect 
ratio may be estimated roughly as the Mach number for which the drag 
coefficient first begins to rise- The rough estimate of the critical 
Mach numbers obtainable by this consideration is not sufficiently accurate 
to warrant comparison of the numerical values with the numerical values’ 
of the critical Mach number obtained in the present paper for thin 
ellipsoids. Comparison of the numerical results is, moreover, not 
warranted inasmuch as the wings of reference 1 did not have an elliptic 
section and furthermore had a rectangular plan form. A qualitative 
comparison may be made, however, between the results of the present 
paper and those of reference 1. The increase in critical Mach number 
with decrease in aspect ratio indicated in figures 1 , 2 , 3 , and 7 of 
the present paper is considered sufficiently large to explain the 
corresponding effect indicated in figure 6 , reference 1 . 

It is mentioned in reference 1 that the Mach number for a significant 
rise in the drag coefficient is approximately 0.1 higher for an aspect 
ratio of 2 than for an infinite aspect ratio. This value is appreciably 
higher than* the increase in critical Mach number due to a decrease in the 
aspect ratio. Since, for low-aBpect -ratio wings, the drag coefficient 
increases only gradually after the critical Mach number is reached, the 
critical Mach number for a wing having low aspect ratio does not indicate 
so critical a change in the flow phenomena as the critical Mabh number 
for a wing having high aspect ratio. It is thought that the smaller rate 
of increase of the drag coefficient for wings having lew aspect ratio is 
due to the fact that, at the critical Mach number, the rate of increase 
with Mach number of the • incremental velocity is less than for high aspect 
ratios, as may be seen from figures 1 to 3 . 
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CONCLUSIONS 


A study by the Prandtl-Glauert method, of compressibility effects 
and. critical Mach number for ellipsoids of various aspect ratios and, 
thickness ratios indicated the following conclusions: 

1. The flow about the unyawed ellipsoid is analogous to that about 
the infinite unyawed cylinder in that sonic velocity is reached simul- 
taneously along a line, that extends across the entire span of the body 
and is normal to the stream direction. 

2. The critical Mach number for a thin ellipsoid may be predicted 
with good accuracy by means of the Prandtl-Glauert method, and the 
accuracy increases with decrease in aspect ratio. 

3« The compressibility effect on the flow about an ellipsoid decreases 
as the aspect ratio decreases. 

4- The three-dimensional relief for ellipsoids is essentially 
independent of the thickness ratio, for thickness ratios from 0-10 to 0 . 20 . 

•• 

■ 5 . Far ellipsoids of thickness ratio 0.20, the critical Mach number 
increases by about 0.04 when the aspect ratio is changed from » to 2j 
for ellipsoids of thickness ratio 0.10 the increase is 0 . 03 . 

6 . The calculated increases in critical Mach number are sufficiently 
large to explain the experimentally observed increases in the Mach number 
at which the drag first begins to rise. 

7* The experimental 1 y indicated reduced rate of drag rise for low- 
aspect-ratio wings at zero lift as compared with that for wings having 
infinite aspect ratio may be explained qualitatively on the basis of the 
results obtained for the three-dimensional relief for ellipsoids. 


Langley Aeronautical laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Ya., April 23, 1946 
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APPENDIX 

THE PELAITDEL-GMIJERT METHOD FOR THREE' -DIMENSIONAL PLOW 


A derivation of the Prandtl-Glauert method for three-dlmpmslonal flow. - 
A “brief derivation of a form of the Prandtl-Glauert method correct for 
three dimensions may he given as follows: A first-order approximation to 

the subsonic compressible flow about a thin body B, the surface of which 
has the equation 


S(s:, y, z) = 0 

may be obtained by finding a solution of the linearized differential 
equation for the potential 9 of the incremental velocities, 

• P 2 Vjoc + 'Pyy + <P 2 z = 0 (Al) 

* 

where the x-axis is in the stream direction and the incremental 

velocities <P X , <P , and 9 are smal 1 compared with the stream 
«y z 

velocity U- At all points on the surface of B, the potential 9 must 
satisfy the boundary condition 

(U + qg S x + 9^ + 9 z S z = 0 (A2) 

which states that the flow is tangential to B- Since B is assumed 
thin, Sj is small - compared with Sy and S z j consequently the second- 

order term 9^5^. may be neglected, and the boundary condition becomes 


VS* + <fyS y + 9 Z S Z = 0 


In order to solve the boundary -value problem given by equations (Al) 
and (A2) in terms of incompressible flow, the following transformation of 
variables is used: 


x 


1 


9' 



(A3) 
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With this transformation equations (Al) and (A2) become, respectively, 


, + cp ' + cp 1 =0 

z 1 v yy Y zz 

■ (Ait) 

+ cp 1 s + cp 1 s = 0 
y y z z 

(A5) 


Equations (Ab) and (A5) are, respectively, the differential equation and 
boundary condition for the potential cp 1 of the incremental velocities 
of an incompressible flow with free-stream velocity' U, in the z ', y, z . 
space, about a thin body B', the surface of which has the equation 


S(0z', y, z) = 0 


The incremental velocities in the compressible flow are thus 
given by 


u = 


= = 72 <P 


. I = T5 U 


V = 


m 1 m' 1 I 

^ m J r 


(A6) 


w = <p 2 = = jw' 


where u, v, and w and u', v', and w' are the incremental 
velocities at corresponding points in the compressible flow about B 
and the incompressible flow about B 1 , respectively. 


The foregoing analysis establishes the Prandtl-Glauert method 
for three-dimensional flow in the following form: The incremental 

velocities at "a point P of a three-dimensional compressible flow 
field about a thin body B may be obtained in three steps: 

(1) The z-coordinates of al 1 points of B are increased by the 
f act or l/p, where 

P = \fl - 1# 

and where the z-azis is in the stream direction. This transformation 
takes B into a stretched body .B 1 . 

(2) The incremental velocities u', v', w', in the direction of 
the z-, y-, and z-azes, respectively, at the point P 1 in the flow 
field of B' corresponding to the point P in the flow field of B 
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are calculated, as though. B ' were in an incompressible flow having the 
same free-stream velocity as the original compressible flow* 

( 3 ) The values u, v, and. w of the incremental velocities at the 
point P in the campressitle flow field, of the orig inal unstretched 
tody are then found, ty the following equations: 


T - i T ' 

Thus far it has teen shown that, through the transformation given 
in this paper, a campressitle flow that satisfies the boundary conditions 
for the tody B is transformed into an incompressible flow satisfying 
the boundary conditions for the stretched tody B 1 * It can be shown 
further that the stream lines of the compressible flow about B are 
transformed into stream lines of the incompressible flow about B 1 * 
Because of this fact, the method has been referred to in the literature 
(for example, reference 4) as the "streamline -analogy" method. 

The proof is obtained simply by applying the transformations (A3) 
and (A 6) to the equations for the streamlines of the incompressible 
flow about the stretched body B 1 


dx 1 dy 1 dz' 

TJ v w 1 


(A7> 


The application of the transformations results in the equation of 
the streamlines for the compressible flow about the body B 


dx _ dy _ dz 
U v w 


(A8) 


Failure for three -dimensional flow problems of the commonly stated 
forms of the Prandtl-Glauert method * - According to the form of the 
Prandtl-Glauert method given by Prandtl (reference 5) and Yon Karm&n 
(reference 6) , the incremental velocities for a compressible flow about 
a thin body B are the same as the incremental velocities of corresponding 
points for incompressible flow having the same free-stream velocity about 
a body obtained by expanding B in the directions normal to the free- 
stream direction by the factor l/p* That is, for bodies of revolution, 
or two-dimensional bodies. 
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According to Gothert 1 s method, however, 

“ (e, M) = \ ~ Oe, 0) (A9) 

u u 

Thus, Prandtl's and Ton Karman's method is valid only if 

0) 

that is, if and only if the incremental velocity far incompressible 
flow about the bodies under consideration is proportional to the thickness 
ratio. This relation is approximately valid for thin two-dimensionai 
bodies, so that the method of Prandtl and Ton Karmen may be expected 
to be valid for two-dimensional flows- The relation is not true in 
general for three-dimensional bodies; for example, for a very thin 
body of revolution the incremental velocity is more nearly proportional 
to the square of the thickness ratio than to the first power. 

Ton Karman approaches the problem by making the transformation 

7 ' = ftr 
z 1 = Pz 
<P' => 

Under this transformation the linearized equation of compressible 
flow goes into Laplace’s equation; however, the transformed boundary 
condition is not satisfied on the surface of the transformed (contracted) 
body but on the surface of- an expanded body. Thus, the boundary condi- 
tion is not satisfied on the boundary but at points near the boundary. 

This procedure is applicable to two-dimensional problems (as, for example, 
in the thin-wing theory, reference lh) , because the velocity increments 
induced by the equivalent line distribution of singularities vary only 
slowly in the neighborhood of the line of singularities. For a body 
of revolution, however, the velocity increments induced by a line of 
singularities go to infinity at the line of singularities; for such 
bodies, accordingly, the location of the point at which the boundary 
conditlon is satisfied is important. 

According to Goldstein and Young (reference 7), "in compressible 
flow the pressure increase at any point of the body is l/(3 times the 
pressure increase in incompressible flow at the same point. " That is. 
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Comparison of this relation with equation (A9) shows that the Goldstein- 
Young method is also valid for two-dimensional problems hut giveB an 
incorrect result for three -dimensional problems. 
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Figure 2- Variation with Mach pumber cf maximum incremental , ve/odfu for P 

ellipsoids having various aspect ratios. Thickness ratio , O.F5. 
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figure 3- Variation with Mach number of maximum incremental velocity for 
ellipsoids having various aspect ratios. Thickness ratio , 0.20. 
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